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A quantum resonance catastrophe for transport through an AC driven impurity 
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We consider the quantum transport in a tight-binding chain with a locally applied potential 
which is oscillating in time. The steady state for such a driven impurity can be calculated exactly 
for any energy and applied potential using the Floquet formalism. The resulting transmission has 
a non-trivial, non-monotonic behavior depending on incoming momentum, driving frequency, and 
the strength of the applied periodic potential. Hence there is an abundance of tuning possibilities, 
which allows to find resonances of total reflection for any choice of incoming momentum and periodic 
potential. Remarkably, this implies that even for an arbitrarily small infinitesimal impurity potential 
it is always possible to find a resonance frequency at which there is a catastrophic breakdown of 
the transmission T = 0. The points of zero transmission are closely related to the phenomenon of 
Fano resonances at dynamically created bound states in the continuum. The results are relevant for 
a variety of one-dimensional systems where local AC driving is possible, such as quantum nanodot 
arrays, ultracold gases in optical lattices, photonic crystals, or molecular electronics. 
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Driven quantum systems appear in many different con¬ 
texts in physics and chemistry At the same time 

there has been remarkable progress in the controlled de¬ 
sign of nanoscale quantum systems with a high degree of 
coherence and tunability. Systems containing just a few 
molecules are promising candidates for the realization of 
electronic components on the sub-silicon scale (molecu¬ 
lar electronics) [lll-[l^. The quantum transfer of par¬ 
ticles between such localized structures can be well de¬ 
scribed by the tight-binding model in order to gain un¬ 
derstanding of the transport mechanisms involved [i3- 
Another versatile realization of near-perfect tight- 
binding models are ultracold gases in optical lattices with 
a great variety of possible geometries [2l|, where tun¬ 
able local impurities periodic driving [Mi , and 

dimensional crossover [23| have also been realized. Tight- 
binding models are also applicable when investigating 
other driven systems such as quantum dot arrays [3|. Fi¬ 
nally, photonic cavities and photonic crystals have been 
used as quantum simulators to fabricate interesting quan¬ 
tum tight-binding systems [HQ. In practical applica¬ 
tions time-dependent effects such as electromagnetic ra¬ 
diation or gate voltages can be used to manipulate the 
transport properties of nanodevices Q . Driven tunneling 
between two quantum wells has been well studied (sj-Q. 
A natural further development is the transport through a 
driven impurity in an extended structure of coupled wells 
with a finite bandwidth. 

In the present Letter we consider the steady state of 
a generic model system, consisting of a one-dimensional 
tight-binding chain for bosons or fermions with a period¬ 
ically varying potential ji at one impurity site (i = 0) 

H = -J y](c|ci+i+ cf+iCi) - ncos{u)t)clco, (1) 

i 

where we have used standard notation and the hopping 


amplitude is denoted by J. This model captures the es¬ 
sential physics of a single one-dimensional band, which is 
useful for the description of corresponding AC driven ex¬ 
perimental setups mentioned above. Higher energy bands 
may give interesting additional effects at very high fre¬ 
quencies, but this will not change the central results in 
this work since the most interesting physics turns out to 
involve excitations close to the band-edge. 

It is well-known that for a static barrier, the phe¬ 
nomenon of tunneling allows transport for any finite po¬ 
tential strength ji. In strong contrast, the time-periodic 
potential considered here turns out to show resonances 
at special driving frequencies where the transmission is 
completely blocked. In fact, we will show that it is always 
possible to find such a finite resonance frequency uj for 
any combination of incoming momentum k and barrier 
strength fi. This implies that even for an infinitesimally 
small periodic perturbation /i there is a complete break¬ 
down of conductance, if the frequency uj is tuned to the 
corresponding resonance. We call this phenomenon the 
quantum resonance catastrophe. 

The goal of this paper is to calculate the transport of 
an incoming particle with a given momentum k and cor¬ 
responding energy e = —2 J cos k, which is the dispersion 
relation of the model in Eq. m away from the impurity. 
Just like in the static case, the transmission coefficient 
can be determined from the steady state solution of the 
Schrodinger equation {H{t) — = 0. Due to 

the periodicity of the Hamiltonian H{t) = H{t 27 t juj) 
it is possible to use the Eloquet formalism 1,113 to 
express any steady state solution in terms of so-called 
Eloquet states |^(t)) = l ^(^))5 where e is the 

quasi-energy of the resulting (d+l)-dimensional eigen¬ 
value equation {H — idt) \ ^{t)) = e |<F(t)) and the Elo¬ 
quet modes |^(t)) = |<F(t + 27r/cc;)) are periodic in time. 
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FIG. 1: Left: Sketch of the model mapped onto a set of static 
coupled chains. The n = 0 chain is locally connected to two 
side-coupled systems of chains with a corresponding chemi¬ 
cal potential ncu. Right: Dispersion relation with the special 
frequencies u = 2J ± e below which the side coupled chains 
n = =f 1 start support unbound solutions for the case e = 0.5 J. 


Since the potential fi only affects a single impurity site 
for all chains, the solutions for j 7 ^ 0 must correspond to 
wave-like states (unbound solutions) for \e-\-nuj\ < 2J and 
bound states otherwise according to Eq. As indicated 
in Fig. [ 1 ] (right) a critical frequency can be defined 

Ldc = 2J + |e|. (7) 

For uj > ujc all chains with n 7 ^ 0 are outside the band 
|e + nuj\ >2J and correspond to bound states. Below the 
frequency uj = uJc the first side coupled chain starts to 
support unbound solutions. A second unbound solution 
starts to appear below uj = 2J — \e\ and so on. 

Let us first consider frequencies uj > uJc with bound 
states for all n 7 ^ 0 of the form 

(8) 


Using the spectral decomposition 

oo 

\m)= E (2) 

n= —oo 

the eigenvalue equation for a Hamiltonian of the form 
H{t) = Ho 2Hi cos{ujt) becomes discrete in the time 
direction 

Ho\^n) + Hi{\^n+l) + l^n-l)) = (e + nc^)|^n). (3) 

A general steady state on the tight-binding lattice is 
given by 

l^n) = ^^0i,ncJ|O). (4) 

j 


where e + nuj = 2Jsign(n) cosh/^^. Inserting these states 
into Eq. dSD, we arrive at a recurrence relation for the 
coefficients Cn for |n| > 0 , 


JnCn = Cn -1 + Cn +1 with 7„ = - 7(e + nui)^ - 4J2, 

/i 

(9) 


where we have defined Co = A cos 0. The solution for this 
second order recurrance relation is fixed up to an overall 
constant by requiring convergence for |n| ^ oo and can 
be solved efficiently numerically. The angle 0 is then 
given by Eq. for n = 0 in terms of those coefficients 


tan^ = 


R 

2uk 




( 10 ) 


The model in Eq. o therefore results in the following 
set of coupled equations 

— + 01,n) ~ ^(0O,n+l + 0O,n-l) = (^ + ^<^)0O,n 

-J(07_i,n+ 0i+i,n) = (e + nu;)0j,n for j 7 ^ 0 (5) 

which effectively corresponds to a static Hamiltonian 
with eigenvalue e for an infinite number of chains labeled 
by n, each with additional overall chemical potential of 
ncj, analogous to a Wannier-Stark ladder [H. The chains 
are coupled to each other only at site j = 0 by a hopping 
term /i /2 as depicted in Fig. [T] (left). Notice that the 
entire problem is symmetric under parity transformation 
j — j, so that solutions are either parity symmetric 
or parity antisymmetric. The parity anti-symmetric so¬ 
lutions obey 00 ,n = 0, Vn, so they do not couple to the 
driving potential and can be ignored. 

Transmission coefficient .— We now want to calculate the 
transmission of an incoming particle with momentum k 
and e = —2 J cosk for the chain n = 0. The parity sym¬ 
metric solution is given by plane waves of the general 
form 


(6) 


where we defined Uk = 2 J sin k as the particle velocity. 

Since the bound states for |n| > 0 do not contribute in 
the transmission, it is now straight-forward to calculate 
the transmission coefficient to be 


T = cos^ 0 = _ - _ 

«2 + (^(C_i-C'i)/2C'o)2- 


( 11 ) 


The transmission obeys T(e) = T(—e). For e = 0 the 
solution becomes symmetric Cn = C-n, which results in 
T(e = 0) = 1 independent of ja for uj > uJc due to Klein 
tunneling [HI. In the following we assume e 7 ^ 0. 

Next we consider lower frequencies uj < uJc^ when un¬ 
bound states also exist for n 7 ^ 0. In this case it is use¬ 
ful to make an ansatz for an incoming wave at energy e 
and transmitted/reflected waves in all unbound channels 
(ioj . In this way it is possible to solve for all parame¬ 
ters. Finally, the total transmission coefficient can again 
be expressed by the solution of the recurrence relation in 
Eq. ®, which now involves the remaining bound states. 
Results .— Using this procedure the exact numerical so¬ 
lution for the transmission coefficient was obtained as 
shown in Fig. [ 2 ] for a given energy e = 0.5 J as a func¬ 
tion of fi and UJ. Perfect transmission T ^ 1 can be 


= Acos(|j|fc - 9). 
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FIG. 2: Exact results for transmission coefficient T as a func¬ 
tion of amplitude /i and frequency uj for an incoming wave 
of energy e = 0.5J. Top: The solid lines (blue) indicate the 
exact values of the resonances (T = 0) while the dashed lines 
are analytical expressions from the first zeros of Bessel func¬ 
tions J±e/uj in Eq. (O at large frequencies and the small /i 
approximations in Eq. d. Bottom: Behavior close to the 
critical frequency cc ~ ccc = 2.5J. Eor uj — ujc the transmis¬ 
sion approaches T (2J — |e|)/8J as /i ^ 0. 


observed for small // or large frequencies. For increasing 
fi there is a sharp drop, however, and at special reso¬ 
nances a vanishing transmission T = 0 can be observed 
(blue solid lines). Interestingly, the transmission then in¬ 
creases again with increasing fi before more resonances 
with T = 0 are reached and so on. This apparent non¬ 
monotonic behavior with potential {i and frequency uj 
can be understood in the high frequency limit, where the 
recurrence relation can be solved analytically. In partic¬ 
ular, note that 7 ^ ^ 2|nci; + e|//i for cj ^ J, so that 
Eq. (|9]) becomes exactly the defining recurrence relation 
for the Bessel functions in this limit. For convergence 
as |n| ^ oc the coefficients therefore can be chosen to be 
Bessel function of the first kind 

Cji ~ (/i/cj). ( 12 ) 

The recurrence relation then approaches Co ~ J±e/uj foi* 
n ^ 0^ from above/below. Accordingly Eq. (pd)|) can be 
approximated ioi uj^J 

^ 2u„ V )■ 

These Bessel functions explain some of the observed fea¬ 


tures of T in Eq. m, namely the oscillating behavior 
with /i and uo and resonances of T = 0 close to the ze¬ 
ros of the Bessel functions J±e/uji which are marked in 
Fig. [ 2 ] as dashed lines (green). Moreover, the behavior 
in Fig. [2] indeed only depends on the ratio /i/cj for large 
frequencies. 

While the description in terms of Bessel functions is 
useful, this does not explain the behavior in the most 
interesting region close to uj ^ ujc where the data shows 
large gradients in the transmission coefficient. In Fig. [2] 
(bottom) the behavior changes dramatically with fre¬ 
quencies just above or below cjc = 2.5J for small driv¬ 
ing potential ji. There is a resonance with T = 0 which 
quickly shifts to smaller /i as the frequency is lowered 
{uj = 2.blJ and 2.5009J) and suddenly disappears com¬ 
pletely once the n = —l chain supports unbound solutions 
{uj = 2A9J). Exactly at the critical frequency uJc = 2.bJ 
the results for ^ 0 show a well-behaved finite value 
T {2J — |e|)/ 8 J, which corresponds to C‘^iICq 7^.2 
and is neither close to unity nor zero. By looking very 
carefully one observes that there is another resonance 
close to yU ~ 11.5J which disappears at ujc- Away from 
these singular points the changes of the transmission are 
very small, however, and appear to be continuous as the 
frequency goes through the critical value. 

Resonances. — It is worth noticing that the resonances 
T = cos^ 0 = 9 ioY UJ > ujc are special points at which the 
coefficient 0o,o = Cq = AcosO vanishes exactly. In this 
case the side coupled systems of the corresponding static 
problem in Eq. ([5]) become decoupled from the chain with 
n = 0 in Fig. [T] (left). Therefore a resonance with T = 0 
for UJ > ujc occurs if and only if the isolated side system 
has an eigenenergy inside the band |e| < 2J. 

This is a remarkable statement since the decoupled side 
chains for n 7 ^ 0 only support bound states outside the 
band. However, due to the local coupling /i between the 
chains one of these energies is pushed inside the band, for 
which T(e) = 0. Such Bound States in the Continuum 
(BIG) were first proposed by von Neumann and Wigner 
for a spatially oscillating potential in the early days of 
quantum mechanics [3l|. Since then BIC’s have received 
exte nsive a ttention in the context of transport phenom- 
The suppression of transmission is closely 
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ena 

related to the Fano-Effect in this case [35, 

In order to illustrate the connection with the Fano ef¬ 
fect the behavior as a function of incoming energy at fixed 
frequency uj = 3 J is shown in Fig. [3] (top). The character¬ 
istic asymmetric lineshape of Fano resonances is clearly 
visible. In analogy to the critical frequency in Eq. 0 , 
it is possible to define a critical energy Cc = cj — 2 J, 
above which the first side chain supports an unbound so¬ 
lution. We observe that a sharp resonance occurs just 
below this energy for small fi which then broadens and 
moves quickly away from this point with increasing fi. 
The lower part of Fig.[3]shows the behavior for larger fre¬ 
quency UJ = 10 J, where the resonances are further apart 
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FIG. 3: Transmission coefficient as a function of incoming par¬ 
ticle energy e at cc = 3J (top) and uj — 10J (bottom). The 
dashed lines in the bottom plot depict the high frequency ap¬ 
proximation in Eq. (O. The inset shows an enlarged region 
of the resonance for /x = 25 J. 


and less pronounced (see inset). In this limit the agree¬ 
ment with the high frequency approximation in Eq. m 
fits reasonably well (dashed lines in Fig. [3]). The drop 
of the transmission T ^ 0 for e ^ 2J occurs due to 
the vanishing of the particle velocity Uk ^ 0 and is not 
connected to any resonance phenomenon. 

As shown in Fig. [3] for fixed ji and uj there is at most 
one energy with T = 0 and in some cases it is not possible 
to find a resonance at all, e.g. for oo = 3J and fi = 10J. 
This is in contrast to the situation of a given energy dis¬ 
cussed above in Fig.[2]where there are always one or more 
resonance frequencies for any value of fi. 

In order to predict the location of the resonances with 
T = 0 we now follow the strategy to consider the eigenen- 
ergies of the decoupled side system in Fig.[T]using Eq. ^ 
with 00,0 = 0- In fact, the two sides for n > 0 and 
n<0 have identical eigenenergies e, due to the symmetry 
transformation (j)j^n Of course there 

are infinitely many eigenenergies, but only those inside 
the band |e| < 2J are of interest. Let us focus on the 
resonances T = 0 close to the critical frequency uj ^ ujc 
in the limit of small ji J. In this case, Eq. m must 
still hold for Co = 0, with ^ 1 and 71 ^ 1. Look¬ 
ing at the first few terms of the recurrence relation it 
becomes clear that the coefficients grow beyond bounds 
unless 1 — 7172 ^ 7i ^ 1- Solving for the frequency at 
which 7172 = 1 we find for the resonance positions 


UJ 


64J([4J-e]2-4J2)’ 


(14) 


which is marked by dashed lines (red) in Fig.[2]and agrees 


well with the exact results. Using the condition 7172 = 1, 
resonances can also be found at fixed uj for small devia¬ 
tions from the critical energy ec, which are again propor¬ 
tional to 

Conclusions .— In summary we have developed a frame¬ 
work to study the transmission across an AC driven bar¬ 
rier connected to leads of finite bandwidth based on the 
Floquet formalism. The resulting transmission coefficient 
T can be calculated exactly and shows versatile tunabil- 
ity with frequency cj, energy e, and impurity strength ji. 
At high frequencies uj ^ J the transmission can be ex¬ 
pressed in terms of Bessel functions, which also appear in 
the description of so-called coherent destruction of tun¬ 
neling. In this limit the oscillations can be averaged to 
form an effective quasi-static hopping Jeff = JJo{ii/uj) 
0, [13 • In our case, a slightly more refined picture 
emerges in terms of Bessel functions with a fractional 
index n = \n\ ± ejuj. 

However, much more interesting effects appear at lower 
frequencies uj ~ uJc = 2J+|e| where there are sharp 
changes in T and a complete breakdown of the transmis¬ 
sion T = 0 may appear even for arbitrarily small barriers 
fi. The explanation of such a quantum resonance catas¬ 
trophe can be found in the dynamically created side cou¬ 
pled chains in Fig. [U which contain bound states for all 
energies outside the band for uj > ujc. The effect of the 
local coupling fi between the chains is to push one energy 
from just above the band into the continuum. Thus effec¬ 
tively a discrete Bound State in the Continuum is formed, 
which is known to have drastic effects on the transmis¬ 
sion [31I based on the Fano effect [HI . There has been an 
interesting proposal to use static side coupled systems to 
engineer Fano resonances [13, but the creation of “vir¬ 
tual” side coupled systems by AC driving considered here 
is even simpler and more versatile than any static design. 
The location of the resulting resonances for /x 0 can 
be predicted rather accurately by Eq. d- It should be 
noted that the width of the resonance also changes dra¬ 
matically near ujc. While this may become the limiting 
factor to observe the quantum resonance catastrophe for 
very small /x, it also presents a unique opportunity for the 
design of switches, where a huge change of transmission 
for small parameter changes near sharp resonances is a 
desirable feature. 

The underlying tight-binding model has long been used 
in condensed matter, but with the recent trend of design¬ 
ing tailored quantum systems, this model has made a re¬ 
vival for the realistic description of corresponding setups 
in molecular electronics, quantum nano dots, and pho¬ 
tonic materials. For ultracold gases in optical lattices 
0 it is now possible to insert localized impurities [22|, 
which play the role of a local barrier that can easily be 
periodically changed using Feshbach resonances. Given 
the large variety of applicable systems it is difficult to 
anticipate which experimental realization is best suited 
to explore the quantum resonance catastrophe predicted 
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